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Abstract We study the recently discovered phenomena of sudden death of entangle- 
ment for a system of two qubits, each of them independently longitudinally damped by 
a reservoir and subjected to a continuous driving. We show that driving produces, in 
the interaction picture, an effective bath that has elements amounting to various extra 
sources of noise (transverse, thermal squeezed, thermal longitudinal). As a result, the 
time of sudden death decreases due to driving, which we verify as well by direct nu- 
merical calculation. We suggest that this phenomenon can be studied systematically 
using superconducting qubits driven by microwave fields. 
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1 Introduction 

Decoherence is an ubiquitous phenomenon in quantum mechanics and any architecture 
of future quantum computers must consider it seriously. Recently, it has been noted 
that already at the level of two qubits, the effect of a weak noise on each qubit is 
not simply additive when it comes to nonlocal properties. This effect, called "sudden 
death of entanglement" [T], means that the concurrence of two qubits has an abrupt 
decay in time rather than the exponential decay typical for single qubits. In this paper 
we consider the case of two qubits coupled independently to two reservoirs, and driven 
coherently and resonantly by two continuous fields. We find that two-qubit states which 
in the absence of the drive would have exponential decay now exhibit sudden death of 
entanglement. 

Although the results in this paper are valid for any type of qubits, we will briefly 
analyze the case of superconducting-qubit architectures, where microwave fields ap- 
plied to each qubit can be used to manipulate the quantum two-qubit state, creating 
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entanglement, CNOT and iSWAP gates, etc. [2]. With the emergence, during the last 
decade, of a plethora of such superconducting devices as serious contenders for the role 
of the quantum bit, the effect of decoherence in these systems should be thoroughly 
investigated. While for the case of single superconducting qubits this has been already 
done to a large extend both experimentally and theoretically, there is yet no system- 
atic experimental study of decoherence in superconducting many-qubit systems. Also, 
circuit boards containing two superconducting qubits could be used as a workbench 
for fundamental studies of decoherence effects. 



2 The model 

We study a general model, similar to that of Ref. [l], which can be adapted in a 
straigthforward way to any species of superconducting qubits. We take ft = 1, and 
we consider two non-interacting qubits which are coupled to two uncorrelated thermal 
reservoirs: 

Hsys = - ^ ^aj, Hres = ^ ^O^jA^A^jA, 
j=1.2 j = l,2 A 

Hs-r = 51 51 (fjA^jA'^^ + 93\aj\cr+^ , (1) 
j=l,2 A 

where ujj' denotes the Larmor frequency of the qubit j, a~ — \ Tj)(ij I i<^j^ = I Ij 

I) is the qubit lowering (raising) operator, and a^^ {(^jx) denotes the bosonic 
creation (annihilation) operator for the reservoir coupled to each qubit, respectively. 
Each qubit is irradiated with an external monochromatic microwave field. The qubit- 
field interaction is taken to be of the form 

Hext{t)= O-j{t)cos{L0jt)G^ , (2) 

with LUj the angular frequency of the field applied to qubit- j, and Qj the qubit-field 
coupling strength (Rabi frequency). 

In the weak external field approximation 3 \\Hext\\ ^ ll-f^saslli which is valid in 
most experiments with superconducting qubits, and zero temperature, the dynamical 
evolution of the qubits' density operator ps is governed by the Born-Markov master 
equation 

ps = -i[Hsys,Ps\ - i[Hext,Ps] + ^Ps]- (3) 

The Liouvillean 

models the longitudinal dampings of the qubits, and Fj is the standard energy relax- 
ation rate. 

The dynamics of the two-qubit system given by Eq. ([3| is rather complex. However, 
it is simplified considerably by transforming into the rotating reference frame, and 
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performing a standard rotating wave approximation (RWA) to eliminate the counter 
rotating terms. The master equation in the rotating frame becomes 

Prf ^-i [Hrf ,Prf]+^ [Prf] , (5) 

with = S^psS, S — exp[i{ijjiaf + tt;2C"2)*/2], and 

+ BS^ 

Hrf = S^Hsys + Hext)S + i-^S 

Here 5j = LUj — lo^ is the detuning of the driving frequency from the corresponding 
qubit transition frequency. 



3 Discussion 

As in other studies on entanglement sudden death in two-qubit systems such as [11[5], 
we assume that the initial state of our system is an "X state" 



/ ao uiQ ' 
6o ^0 
zo CO 
do. 



(7) 



However, despite the fact that the "X state" is a rather general form of a density matrix, 
due to the existence of driving fields, the evolution under Eq. ([5]) does not preserve 
this form, resulting in a rather complicated kinetic equation for the density matrix 
elements. Thus, we need to find a way to reduce the number of coupled differential 
equations. 

For the sake of simplicity, we assume on-resonance driving [5j = 0) and time 
independent Qj (rectangular pulses), then we move to the interaction picture by means 
of p = exp{iHrft)prf exp{—iHrft). The equation of motion for p is thus 

i=i,2 

where 



o-j = e ' ^ ' a - e ' ' ' = -- sm{njt)(jj + sm I -|- 1 aj + cos I I o-p,) 

~+ ifijCr°'t/2 + -il7,CTjt/2 ^ ■ i n ,\ z , 2 ( + , . 2 ( -r t n\ 

cr J = e ' ' ' aJe ^ ^ ' = - sm(i2jt)crj + cos I I °" J + I I cTjiW) 

Substituting ^ and (|10p into ([8]), the master equation in the interaction picture 
is explicitly given by (for notation simplicity, we drop the tilde over p from now on) 
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+ rj[3 + Acos{njt)+cos{2njt)](2a-pa+ -a+a-p- pa+aj) 

+ — ^ rj[3' AcosiOjt) +cos{2njt)][2a+paJ - aja+p- paja^) 
j=l,2 

+ cos(2r2jt)] + cj-paj + pa] - p) . (11) 

It is easy to verify that when the external fields are very weak, Qj <C Fj, even in 
a relatively long time scale r oc max[/y"^], Eq. (|lip is reduced to 

A more realistic situation is the so-called secular limit Oj ^ Fj [6j. For instance, 
in experiments with quantronium circuits |7], the Rabi frequency fij is usually larger 
than 100 MHz, whereas the relaxation rate Fj < 2 MHz. In this case, we may neglect 
all the oscillating terms with frequencies Oj and 2J7j in Eq. (|ll|l . and arrive at a simple 
equation of motion 

P~ ^ ((7+pa+ + a^paj^ + a] pa] - p) 

J = 1.2 

+ E T# {'^"Jp^t + 2'^/^^7 - 4^7^ - ^74^^ - p'^t^J - p^J^t) ■ (13) 

J = l,2 

We note that Eq. (|13|l can be regarded as describing losses caused by a combination 
of three uncorrelated effective noise sources of each qubit: 

transverse thermal reservoir 

^aj[p\ = ^{'T]pa]-p); (14) 
longitudinal high temperature thermal reservoir 

^bj[p] = + 1) (2'^7^'^/ ^ '^t'^-fP ^ P'^t'^j) 

with fij + 1 ~ fij and 'yjnj/2 ~ Fj/8; 
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longitudinal squeezed vacuum reservoir [8] 

^cjIp] = Y [cosh^(rj) (2(7^" pa + - a+a^p-pa^a^) 

+ sinh^(rj) (2cr+paj7 - cr^^cr+p - pcrj^cr+)] 
— Atj sinh(rj) cosh(rj) 

X (^e-'"' a- pa J + e^'^ <T+pa+) , (16) 

with 9j — ±7r being the reference phase for the effective squeezed field, rj being 
the squeeze parameter which fulfills sinh(rj) « cosh(rj) (large squeezing limit) and 
Hj sinh(rj) cosh(rj) ~ -Tj/S. 
In other words, 

p= + [-"i + 

3 = 1,2 

Therefore even without further calculations, we can predict that the external driv- 
ing fields will enhance the sudden death effect, i.e., it will destroy the long-lived entan- 
gled states investigated by Yu and Eberly [1], due to the raising of reservoir's effective 
temperature and combination of noises [9]. 



4 Sudden death time 

Because of its invariance under local unitary transforms, the concurrence of the system 
can be estimated by solving Eq. p3|l , which is much simpler than solving ([SJ since the 
X states remains the X form in the evolution given by (|13|l . Since p(0) — Ps(0), we may 
set 

(a{t) w{t)\ 
bit) z(t) 
z*it) cit) ' ^^^> 
w*{t) d{t) J 

and by substituting it into (|13p we obtain the following kinetic equation for those 
non-zero density matrix elements: 



a = 


3[A(c- 


-a)+r2{b- 


a)]/8. 


b = 


3[ri(d- 


-b) + r2(a- 


b)]/8, 


c = 


3 [A (a 


-c) + r2{d- 


c)]/8, 


d = 


3[A(fc- 


-d) + r2{c- 


d)]/8. 


z — 


[A(»* 


~5z) + r2{w 


"52)]/8, 


w - 


-- [A(.* 


-5w) + r2{z 


-5w)]/8. 



(19) 

By defining rjj = exp(— J^jt/4), the solution of Eq. p9|) has the form: 



a = 
b = 
c = 



[ao (i + r?i + ryl + mm) + bo {i + vi - V2 - mm) + CO (i - ??i + 7?2 - mm) +do{i~vi -V2 + V1V2)] /4, 

[ao {I + V1-V2- vhl) +bo{l + Vl+V2+ vhl) +c^{l-Vl-V2+ vM) + do (l - r?i + 772 ~ V1V2)] /4, 
[ao (1 - ryi + r?2 - mV2) + bo {l - Vi - 1)2 + r?!??!) + co (l + m + V2 + V1V2) + do {l + Vl - V2 - mV2)] /4, 
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d = [ao (l - ,,3 _ ^3 _^ ^3^3^ + bo {l - Vl + m - mV2) + Co {l + Vl - V2 - Vim) +do{l+Vl+V2+ vM)] /4, 

z = ?7i% ["'o ii + vi-V2- mv2) + zoii + 111 + m + mm) + ^o (i - -ni + m - mm) + zo{i-m-m + mv2)] /4, 

w = V1V2 [wq (1 + VI +V2 + V1V2) + zo (1 + ?7i - ?72 - V1V2) + mo (1 - '71 - % + V1V2) + zq {1 - Vl + V2 - V1V2)] /4. 

(20) 




Fig. 1 F and G as functions of rj and a. 



To obtain the sudden death time, we follow Sec. Ill in [S]. For even greater sim- 
plicity, we assume Fi = r2 = F, thus rji = ri2 = v, and < ry < 1. The two equations 
for the sudden death time, 

F{t) = \z{t)f - a{t)d{t) = 
and G(t) = h(t)P - b{t)c[t) = 0, (21) 

are polynomial equations of degree 12 of r], which can not be solved straightforward. 
However, as illustrated in Fig. [T] by plotting F and G vs -q for an initial Yu-Eberly 
state 

/l~a 0^ 



110 
110 
V a, 



(22) 



we find that both functions are continuous in the range < < 1, G is always negative 
in this range, and F becomes zero at 77 « 0.9 for all values of a, which indicates the 
sudden death time tgfi ~ 0.42/F for all a. As we discussed in the preceding section, the 
external microwave fields enhance the disentanglement of long-lived entangled states 
with a < 1/3. 



5 Numerical verification 

In Sec. [S] to derive the analytically solvable master equation (|13p . we have made the 
approximation of neglecting all the oscillating terms. Now, we numerically check the 
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n = 



n = o.5r 




a = 5r 



a = 5or 





Fig. 2 (a) Time evolution of concurrence for the Werner state in the rotating frame (with 
four different driving amplitudes) obtained by solving numerically the initial master equation 
Eq.iO. One can notice that at low Rabi frequencies [Q = 0, 1? = 0.5r) the sudden death 
effect appears only for some of the states, as expected from [l] and Eq. II12I I. However, in the 
regime of larger Rabi frequencies {Q = 5r, Q = SO-T), the sudden death effect is enhanced, 
(b) Time evolution of concurrence for the Werner state in the interaction picture and the limit 
of large Rabi frequencies (i7 = SO-T) , obtained by using the analytical solution of the master 
equation Eq. ((13}, namely Eq. ll20t and Eq. ((25J. 



validity of this approximation. We consider that initially the two qubits are in a Werner 
state g] 



Pw 



-/ 

II (l + 2/)/2 (l-4/)/2 
3 I (l-4/)/2 (l + 2/)/2 

i-f, 



(23) 



8 



with the fidehty 0.25 < f < 1. 

We first calculate the concurrence in the rotating frame with the original formula 

C(p^/) =max{0,Ai-A2-A3-A4}, (24) 

where AjS are the eigenvalues of ^ ^JTrJ Pr j ^JTrJ in decreasing order, with p^j = 
(a^ ® a^)f3*j:{a^ ® cr^). Fig. [2ja) shows the time evolution of C{prf) obtained by 
numerically solving the master equation ([SJ . For simplicity, we take the Rabi frequencies 
f2i = ^2 = ^, and still keep -Ti = ^2 = r. In the interaction picture, p remains in the 
X form during its time evolution, therefore a simplified formula [4l[5] 

C{p) = 2max{0, \z{t)\ - y/a{t)d{t), \w{t)\ - y/b{t)c{t)} (25) 

is applied for computing the concurrence. Fig. [2jb) illustrates the time evolution of 
concurrence for the Werner state obtained by direct substituting (|20p into (|25p . 

By comparing Fig. [2jb) with the lower subplots of Fig. [2ja), one can find that 
the approximation of neglecting oscillating terms is very good for large O. From the 
upper subplots of Fig. [2][a) one can also find that the enhancement of sudden death 
can happen even when O < F. 

6 Conclusion 

In this paper we have studied the entanglement of two non-interacting or interacting 
qubits driven by external on-resonance microwave fields, in the presence of uncorrelated 
zero temperature longitudinal noise sources. We have found that in the secular limit, 
the behavior of the concurrence is very similar to that at finite temperature and under 
both longitudinal and transverse dampings. Although the case of two qubits coupled 
to transverse noises was not considered in this paper, they can be easily investigated 
in the same way as for longitudinal noises. In the secular limit and in the interaction 
picture, the transverse noises will contribute with the following Liouvillean 

3 = 1-2 

-a+ pa+ - aj paj) . (26) 
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